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ABSTRACT 

We calculate the chiral odd quark distributions and the corresponding structure func- 
tions Iit{x,Q 2 ) and hi(x,Q 2 ) within the Nambu-Jona-Lasinio chiral soliton model for 
the nucleon. The Q 2 evolution of the twist-2 contributions is performed according to the 
standard GLAP formalism while the twist-three piece, Iil(x), is evolved according to the 
large Nc scheme. We carry out a comparison between the chiral odd structure functions 
of the proton and the neutron. At the low model scale (Qq) we find that the leading 
twist effective quark distributions, f[ q \x, Qq), g[ 9 \x,Ql) and h^\x,Ql) satisfy Soffer's 
inequality for both quark flavors q — u,d. 
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1 Introduction 

There have been a number of recent investigations into the chiral odd structure func- 
tions of the nucleon. As in the case of the polarized structure functions there are two 
quantities of interest at leading twist: The transverse spin chiral odd structure function 
hr{x,Q 2 ) and the longitudinal spin chiral odd structure function }il(x,Q 2 ). Within the 
context of the operator product expansion (OPE) the analysis in terms of twist reveals 
that the transverse chiral odd structure function hr{x, Q 2 ) is purely twist-2, while the 
longitudinal structure function Hl{x, Q 2 ) contains both twist-2 and twist-3 contributions. 
Accordingly, the decomposition of h^x^Q 2 ) into twist-2 and twist-3 (h L (x,Q 2 )) pieces 
is given by 



As a reminder we note that the kinematics are defined such that q denotes the momentum 
transferred to a nucleon of momentum p. In the Bjorken limit, i.e. Q 2 = —q 2 — ► oo with 
x = Q 2 /2p -q fixed, the leading twist contributions to the nucleon structure functions 
dominate the 1/Q 2 expansion. The additional and important logarithmic dependence on 
Q 2 , which is associated with soft gluon emission, is included via the evolution program of 
perturbative quantum-chromo-dynamics (QCD). 

While the chiral odd structure functions are not directly accessible in deep inelastic 
lepton nucleon scattering (DIS) there is the well known proposal at RHIC to extract the 
quark transversality distributions (x, Q 2 ) (a being the flavor index) from Drell-Yan 
dilepton-production resulting from transversely polarized proton beams Unfortu- 
nately dilepton production processes are difficult to extract from proton-proton collisions 
as the purely hadronic processes dominate. Furthermore this experiment will provide only 
the product of the chiral odd distributions for quarks and antiquarks. As the latter are 
presumably small these flavor distributions are not easily measurable in the Drell-Yan 
process. In the light of these disadvantages it has recently been pointed out that the 
transversality distributions may also be measured in the fragmentation region of DIS [0 . 
The key observation is that these distribution functions can be extracted from an asym- 
metry in the two meson production in the special case that this two meson state (like 
7r + 7r~) is a superposition of different C-parity states, as e.g. a and p. Then the phases in 
the final state interactions do not vanish on the average and the differential cross section 




(1) 
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is proportional to the product of chiral odd distributions and the interference fragmenta- 
tion functions. The latter describe the emission and subsequent absorption of a two pion 
intermediate state from quarks of different helicity. In case these fragmentation functions 
are not anomalously small the chiral odd distribution functions can then be obtained 
from DIS processes[] like eN — > e'n + 7i^X with the nucleon N being transversely polar- 
ized. Assuming isospin covariance for the fragmentation functions these DIS processes 
will provide access to the charge squared weighted chiral odd distribution functions 0. 
Such processes should be measurable in the transversely polarized target experiments at 
HERMES. Knowledge of the chiral odd structure functions will serve to complete our 
picture of the spin structure of the nucleon as they correspond to the distribution of the 
quark transverse spin in a nucleon which is transversely polarized ||. With these data 
being expected in the near future it is, of course, interesting to understand the structure of 
the nucleon from the theoretical point of view. As we are still lacking a bound state wave 
function for nucleon in terms of quarks and gluons, i.e. computed from first principles in 
QCD, it is both mandatory and fruitful to investigate these chiral odd flavor distributions 
and their charge weighted average nucleon structure functions within hadronic models of 
the nucleon |, §, § % g g 0. 

In the context of the spin structure of the nucleon chiral soliton models are particularly 
interesting as they provide an explanation for the small magnitude of the quark spin 
contribution to the proton spin, i.e. the vanishingly small matrix element of the singlet 
axial current [11]]. In these models the nucleon is described as a non-perturbative field 
configuration in some non-linear effective meson theory [|l^, [L3| Unfortunately in 
many of these soliton models the evaluation of structure functions is infeasible due to the 
highly non-linear structure of the current operators and the inclusion of higher derivative 
operators which complicates the current commutation relations. However, it has recently 
been recognized that the soliton solution |14j] which emerges after bosonization [[ilj of 
the Nambu-Jona-Lasinio (NJL) |16[ chiral quark model can be employed to compute 
nucleon structure functions [|17], [18| . In order to project this soliton configuration onto 
nucleon states with good spin and flavor a cranking procedure must be employed [|1^, [19 



which implements significant 1/Nq contributions (Nq is the number of color degrees of 
freedom.). When extracting the structure functions from the NJL chiral soliton model 
the full calculation which also includes effects of the vacuum polarized by the background 



lr The relevant fragmentation and distribution functions depend on different kinematical variables: the 
two meson state momentum fraction and the Bjorken variable, respectively. 
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soliton is quite laborious. In addition we are still lacking a regularization prescription 
of the vacuum contribution to the structure functions which is derived from the action 
functional and which yields algebraic expressions for their moments which are consistent 
with those for the static nucleon properties. Fortunately it is known that the dominant 
contribution to static nucleon properties stems from the single quark level which has the 



lowest energy eigenvalue (in magnitude) and is strongly bound by the soliton |L4| . This is 
particularly the case for spin related quantities. Hence it is a reasonable approximation to 
consider only the contribution of this level to the structure functions. In the proceeding 
section the NJL chiral soliton model together with the above mentioned approximation, 
which we will call valence quark approximation^ will be described in more detail. 

The NJL model for the quark flavor dynamics incorporates spontaneous breaking of 
chiral symmetry in a dynamic fashion. Hence the quark fields which built up the soli- 
ton self-consistently |2(| are constituent quarks with a constituent quark mass of several 
hundred MeV. Keeping this in mind we calculate both the effective constituent quark dis- 
tributions and in turn the corresponding leading twist contributions to nucleon structure 
functions (cf. eq (§)) at a low scale Qq. In the language of Feynman diagrams the DIS 
processes are described by a constituent quark of the nucleon absorbing a quanta of the 
external source. In the Bjorken limit the quark then propagates highly off-shell before 
emitting a quanta of the external source. The intermediate quark may propagate for- 
ward and backward. Hence the complete structure functions acquire contributions from 
both distributions where the intermediate constituent quark moves forward and back- 
ward. We will focus on nucleon structure functions which are defined as the sum over the 
charge-weighted flavor distributions || 

^ a 

in analogy to those of the chiral even spin polarized and unpolarized nucleon structure 
functions || [§. Here a represents a quark label, while (±) refers to the forward (+) 
and backward (— ) propagating intermediate constituent quarks. Furthermore e a denotes 
the charge fraction of the considered quark flavor a. The complete chiral odd structure 
functions are finally obtained as the sum 

h T /L{x, Ql) = h£/ L {x, Ql) + h ( '/ L {x, Ql) . (3) 



2 This notation refers to the valence quark in the NJL chiral soliton model and should not be confused 
with the valence quark in the parton model. 
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The calculation of the flavor distributions hju L in the valence approximation to the NJL 



chiral soliton model [[17], is summarized in section 3. 

Further it is important to note that when considering model structure functions the 
OPE implies that the initial conditions, fi 2 = Q 2 ,, for the evolution is, a priori, a free 
parameter in any baryon model [^TJ. For the model under consideration it has previously 
been determined to Qq « 0.4GeV 2 by studying the evolution dependence of the model 
prediction for the unpolarized structure functions |I7|]. In a subsequent step to compute 
the chiral odd structure functions we employ a leading order evolution program |], [| to 
obtain the chiral odd structure functions at a larger scale, e.g. Q 2 ~ 4GeV 2 relevant to 
the experimental conditions. This evolution program incorporates the leading logarithmic 
corrections to the leading twist pieces. The evolution procedure as applied to our model 
structure functions will be explained in section 4. 

The numerical results for the chiral odd structure functions are presented in section 
5 while concluding remarks are contained in section 6. Technical details on the model 
calculations and the QCD evolution procedure are relegated to appendices. Let us also 
mention that there has been a previous calculation of hx{x,Ql) [|1(| which, however, ig- 
nored both the projection onto good nucleon states and the QCD evolution. Furthermore 
in that calculation an (arbitrary) meson profile was employed rather than a self-consistent 
soliton solution to the static equations of motion. 



2 The NJL-Model Chiral Soliton 

Before continuing with the discussion of the chiral odd structure functions, we will review 
the issue of the chiral soliton in the NJL model. 



The Lagrangian of the NJL model in terms of quark degrees of freedom reads [T5[ 

C = qffl - m°)q + 2G NJL £ ([q T -q) 2 + (g^ 7 5<?) 2 ) . (4) 

Here q, rh° and Gnjl denote the quark field, the current quark mass and a dimension- 
ful coupling constant, respectively. This model is motivated as follows: Integrating out 
the gluon fields from QCD yields a current-current interaction mediated by one gluon 
exchange to leading order in powers of the quark current. Replacing the gluon mediat- 
ing propagator with a local contact interaction and performing the appropriate Fierz- 
transformations yields the Lagrangian (||) in leading order of 1/Nc wnere Nc 



refers to the number of color degrees of freedom. Although only a subset of possible non- 
perturbative gluonic modes are contained in the contact interaction term in eq (f|) it is 
important to stress that gluonic effects are contained in the model (f|). Furthermore the 
NJL model embodies the approximate chiral symmetry of QCD and has to be understood 
as an effective (non-renormalizable) theory of the low-energy quark flavor dynamics. 

Application of functional bosonization techniques |15j to the Lagrangian (^) yields the 
mesonic action 



A = Tr A log(D) + — — / d*x tr (m° (m + M f ) - MM ] ) , (5) 
D = iQ - (M + M f ) - 75 (m - M f ) , (6) 



where M = S + iP comprises composite scalar (S) and pseudoscalar (P) meson fields 
which appear as quark-ant iquark bound states. For regularization, which is indicated by 



the cut-off A, we will adopt the proper-time scheme |p4j| . The free parameters of the 
model are the current quark mass m°, the coupling constant Gnjl and the cut-off A. 
The equation of motion for the scalar field S may be considered as the gap-equation for 
the order parameter (qq) of chiral symmetry breaking. This equation relates the vacuum 
expectation value (M) = ml to the model parameters m , Gnjl and A. For apparent 
reasons m is called the constituent quark mass. The occurrence of this vacuum expectation 
value reflects the spontaneous breaking of chiral symmetry and causes the pseudoscalar 
fields to emerge as (would-be) Goldstone bosons. Expanding A to quadratic order in P 
(around (M)) these parameters are related to physical quantities; that is, the pion mass, 
m n = 135MeV and the pion decay constant, f n = 93MeV. This leaves one undetermined 
parameter which we choose to be the constituent quark mass ||15|| . 



The NJL model chiral soliton [14, 20 1 is given by a non-perturbative meson configu- 



ration which is assumed of the hedgehog type 

Mu(x) = m exp (it ■ xQ(r)) . (7) 

In order to compute the functional trace in eq (||) for this static configuration we express 
the Dirac operator (|6|) in terms of a Hamiltonian operator h, i.e. D = if3(d t — h) with 

h = a. ■ p + m (3 exp (ry 5 r • xQ(r)) . (8) 

We denote the eigenvalues and eigenf unctions of h by e M and respectively. Explicit 
expressions for these wave-functions are displayed in appendix B of ref |14}|. In the 
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proper-time regularization scheme the energy functional of the NJL model is found to be 



19, 14 



rnr^T N C / \ \ N C f°° d,S ^ / 2 \ 

2 2 J1/A2 V47TS 3 V v 7 

+ m 2 J 2 J d 3 r (1 - cosQ(r)) . (9) 

The subscript "v" denotes the valence quark level. This state is the distinct level bound 
in the soliton background, i.e. —m < e v < m. The chiral angle, O(r), is obtained by 
self-consistently extremizing E[Q] [pCf . 



States possessing nucleon quantum numbers of spin and isospin are generated by ele- 
vating the rotational zero modes to time dependent large amplitude rotational fluctuations 
about the hedgehog field [T3| 



M(x,t) = A(t)M n (x)A\t) , (10) 

which introduces the collective coordinates A(t) G SU{2). Substituting the ansatz fllCf ) 
into the action functional (|5|) and expanding jB5] in the angular velocities 

2A\t)A(t) = ir-Q (11) 

to quadratic order yields the Lagrange function for the collective coordinates. Upon 
canonical quantization the angular velocity f2 is substituted by the nucleon spin operator 
J = a 2 0, with a 2 being the moment of inertia [ |H | . The eigenf unctions of the resulting 
Hamiltonian are the Wigner D-functions 

<^> = ^ D% £-jM) . (12) 

with J 3 and J 3 being respectively the isospin and spin projection quantum numbers of 
the nucleon. The nucleon matrix elements of the collective rotations are obtained via 
(N\\x{TiATjA<)\N) = -(8/3)(N\IiJj\N) [0. This approach to generate nucleon states 
from the hedgehog corresponds to the cranking technique in nuclear physics ||25|| . 

Expectation values of bilocal quark-bilinears appearing in the evaluation of nucleon 
structure functions are expressed as (regularized) sums over bilocal and bilinear combi- 
nations of all eigenfunctions *Sf ^ including the Dirac sea states. In practice this is quite a 
painful task, in particular when cranking corrections (|T0|) are included. Also the problem 
of regularization is not consistently solved. Fortunately it turns out that the dominant 
contributions (> 80%) to static nucleon properties (which are moments of the structure 
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functions) stems from the distinct valence level ^> v fl4| . It is therefore reasonable to ap- 
proximate the relevant bilinears by their valence quark contribution. In order to obtain 
quark distributions of the nucleon and the corresponding nucleon structure functions (c.f. 
0), rather than soliton structure functions the cranking contribution to the wave-function, 
which is induced by the collective rotation A(t), must be included. That is, the valence 
quark wave-function employed to approximate the bilinears in the structure functions 
reads 

* v (jM) = e-^A(t)U v (x) + W ^(^) (/i|T l° |v) } = : e-^A^ix). (13) 

Here ip v (x) refers to the spatial part of the body-fixed valence quark wave-function 
with the rotational corrections included and ^ = (ji\x) are eigenfunctions of the Dirac 
Hamiltonian (||). This replacement of the bilocal and bilinear quark fields when computing 
nucleon structure functions defines the valence quark approximation. 

3 Chiral Odd Structure Functions, hr(x) and hi(x) in 
the NJL model 

Here we present the major topic of this paper, namely the calculation of the twist-2 
and twist-3 chiral odd structure functions in the NJL chiral soliton model. Like their 
deep inelastic chiral even (un)polarized counterparts, the chiral odd structure functions 
are computed as Fourier transformations of nucleon matrix elements of bilocal quark 
operators on the light-cone The key features of the relevant light-cone kinematics 
are given in Appendix A. 

We begin by listing the forward propagating intermediate quark (+) contribution to 
the chiral odd nucleon structure functions. Before, however, straightforwardly transcrib- 
ing the expressions from appendix A we must recall that the soliton represents a localized 
field configuration. Therefore a collective coordinate x Q is introduced which parameter- 
izes the position of the soliton. This collective coordinate is employed to generate states 
with good linear momentum |[2G|| . When computing matrix elements between states of 
identical momenta one is essentially left with an integration over x . In the nucleon rest- 
frame (RF) the contribution of the forward moving intermediate quark to the chiral odd 
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structure functions may therefore be expressed as{] 



S7T J 



Mi, 

7f' 



I d 3 x (S ± \^ + (C-x ) 1±lB Q 2 ^ + {-x )\S ± )^__ 



(14) 



For convenience we have omitted the subscript "v" for the valence quark wave function. 
Note that £ refers to a four vector which in light-cone coordinates reads (£ + , £~, This 
coordinate enters the light-cone variables via ^ = (t ± z)/ y/2. Also, the notation S± is 
synonymous for the spin being perpendicular to the coordinate z. On the other hand for 
the longitudinal counterpart 

2MV2 



h[ +) (x) = N c —^- /drexpHC 

l07T J 



Mx, 
72' 



J d 3 x (S z \^-x ) lol ,Q 2 ^_(-x ) 



- Zo)7o75Q 2 ^+(-*o)|S 2 > 5+=€±= o (15) 



the spin is aligned with the z-axis. The "good" and "bad" light-cone components of 
the quark wave functions are the projections *ff± = P±9, with P± = |7 =F 7 ± being the 
corresponding projections operators. Above, Q = diag(2/3, — 1/3) refers to the matrix 
containing the quark charge fractions and the zero momentum nucleon states are given 
by, \p = 0, S) = [(27r) 3 2M] 5 \S). Introducing Fourier transforms for the spatial part of 
the valence quark wave functions (cf. eq (|13"D ), 



yields, 



h 



(+), 



-a 

>/2 J 



M 



c 



V2 



TV" 



d 2 p±dp 3 
2vr 2 



exp 



. V2 



df; p 2 dp d(cos9) d<p exp 



(Mx 



pcosf 



(16) 



V2 J 

x(s ± \ijUph±^Q 2 Mp)\s±) , (17) 



and 



N, 



M 



c 



y/2- 



71" 



d£ p 2 dp d(cos( 



x(s z \ipl( P ) lol5 g 2 i>4p) - ^-(p)ioi 5 q:'Mp)\Sz) • (18) 

for nucleon 



exp 



— i£ (Mx — e v + pcos i 



3 The following expressions constitute a generalization of Jaffe's original definition 
structure functions. 
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In anticipation of the decomposition ([13]) the square of the charge operator is redefined 
to 



,2 5 1 



Q =Y8 1+ Q Dsm - (19) 
Here — ~ tr (riA{t)TjA^ denotes the adjoint representation of the collective rotation 
which is defined in eq (PH|). The integrals over £~ and 9 enforce both the constraint 
cos 6 = (e v — Mx)/p and the lower bound, p m \ n = \Mx — e v | on the p integration. This 
results in the forward moving quark contributions to the transverse and longitudinal chiral 

odd nucleon structure functions 

2M 



/4 +) (x) = N c / pdp d<j>{S ± $l(p)w 5 Q?j> + (p)\S ± ) cose= ^ (20) 

7T J Pmin P 



and 



= N c — r pdpd^SMiPholsQ^-iP) 

-^-(pholsQ 2 Mp)\S z ) cose =^ ■ (21) 

v 

In order to obtain the full structure functions the contribution of backward moving quarks 
hj> L must be considered as well. These contributions are easily obtained from h^ L (—x) 
by reversing the appropriate signs in equations (|H]) and (|TEp. The two contributions may 
be comprised inf] 

, , - A/T roc 

hg\x) = ±N C - P d P d<t>(S ± \^(p )(lTa 3 ) 1±l5 Q 2 iP(p T )\S ± ) I ^ (22) 

r rain 

and 

hf\x) = ±Nc-J^ pdpd^(S z \^(p T )a 3l0l5 Q 2 ij(p T )\S z )\ cos ^ (23) 

where p^ in = \Mx±e v \ and cos^ = (Mx±e v )/p and ip(p±) = ip(p, cosO^, <p). Finally we 
summarize our results by decomposing the proton structure functions into their (iso) scalar 
and vector components, 

h T ( x ) = /4= + °(x) + h^Kx) + (h^{x) + 4 = i(x)) , 

h L (x) = hfflx) + h{^(x) + (h I L =°(x) + h I =l(x)) . (24) 

The isoscalar piece (I = 0) originates from the unit matrix in the decomposition flT§|) 
while the isovector part (/ = 1) stems from the terms involving the collective coordinates. 
The explicit expressions for the structure functions (24) in terms of the static quark wave 
functions are computed in Appendix B. 



4 We have used that the valence quark level has positive parity, i.e. under p — » p we find ijj(-p) = 
7V(p). 
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4 Projection And Evolution 



We consider that our model approximates QCD at a low scale Qq. In oder to compare the 
predicted structure functions with data they must be evolved to a (larger) Q 2 commen- 
surate with experimental conditions. A direct comparison with data gathered at a low 
scale cannot be made as the latter structure functions contain sizable contributions from 
higher twist. Thus we evolve the chiral odd model structure functions of the preceding 
section utilizing the results of perturbative QCD. 

In the soliton approach the baryon states are built from localized field configurations. 
In fact, these states do not carry good four-momentum. Therefore the calculated structure 
functions (cf. Figs. [I], 3 an d do not vanish exactly for x > 1 although the contributions 
for x > 1 are very small. 




-0.25 



-0.25 



Figure 1: The valence quark approximation of the transverse chiral-odd nucleon distri- 
bution function as a function of Bjorken-x for the up and down quark flavor content in 



the rest frame. For comparison also the model calculation [y| for the twist-2 polarized 



structure function gi{x,Ql) is shown for the respective flavor channels. Two values of 



the constituent quark mass are considered: m 
(right panel). 



400MeV (left panel) and m = 450MeV 



The calculation of nucleon structure functions in the Bjorken limit, however, singles 
out the null plane, £ + = 0. This condition can be satisfied upon transformation to the 
infinite momentum frame (IMF) even for models where the nucleon emerges as a (static) 
localized object p8[. For the quark soliton model under consideration this transformation 
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corresponds to performing a boost in the space of the collective coordinate Xq, cf. eq (|14|). 



Upon this boost to the IMF we have observed |29| that the common problem of improper 
support for the structure functions, i.e. non-vanishing structure functions for x > 1, is 



cured along the line suggested by Jaffe |30[ some time ago. The reason simply is that the 

Lorentz contraction associated with the boost to the IMF maps the infinite line exactly 

onto the interval x £ [0, 1[. In addition we have observed that this Lorentz contraction 

effects the structure functions also at small and moderate x. Incorporating these results 

for the general set of leading twist structure functions within the NJL-chiral soliton model 

yields the following form for the forward and backward moving intermediate quark state 

contributions to the chiral odd transverse spin structure function, (x,Q 2 ), 

(±) M f°° 
h> T >(x) = ±N C — -/ pdpdip 

TT[l — X) Jp min 

x(N\ft(p T )(lT« 3 )7±7 5 Q 2 ^( P T )\N) M H i- X)±ev ■ (25) 

cosf/=— 1 

V 

In general the resulting relation between structure functions in the IMF and the rest frame 
(RF) reads 

/imf(x) = 9 [ 1 _~/ ) /rf( - ln(l - *)) • (26) 
Of course, in the context of the chiral odd structure functions /rf is to be identified with 



the expressions in eqs ([22|,[2l||4]) . As will be recognized shortly the solution to the proper 
support problem is essential in order to apply the evolution program of perturbative QCD. 
The chiral odd and polarized structure functions resulting from this transformation are 
shown in figure 

In order to include the logarithmic corrections to the twist-2 pieces of the chiral odd 
structure functions we apply the well-established GLAP procedure For the transverse 
component Ht{x, Q 2 ) this is straightforward as it is pure twist-2. For the longitudinal 
piece Hl{x, Q 2 ) one first has to extract the twist-2 component through Ht{x, Q 2 ) namely, 
h ( i\x,Q 2 ) = 2x Si dyh T (y,Q 2 )/y 2 . 

We simultaneously denote by the twist-2 parts of hx and h L . To leading order (in 
a QCD{Q 2 )) the variations of the structure functions from a change St of the momentum 
scale is given by 

(x, t + St) = ( Xjt ) + rf/ * (2) , (M) St , (27) 

at 

where t = log {Q 2 / -^qcd) ■ The variation fl2~7|) is essentially due to the emission and 
absorption of soft gluons. The explicit expression for the evolution differential equation 
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Figure 2: Same as figure [I] in the IMF fl26|). 



is given by the convolution integral, 

dh (2 \x,t) _ a QCD {t) 



dt 



2vr 



Cr{F) 



dy 

y 



P q \ (y) h (2) 



X 



;28) 



where the leading order splitting function |32| g is given by, 

2 



2 + -6(z-l) 



(29) 



and Cn(f) = (rvj — /2nf for n./ active flavors, aQCD(t) = 47r/ [b log (Q 2 /A 2 )] and 
&o = (ll-Afc — 2rif)/3. Employing the "+" prescription yields for three light flavors and 
N c = 3 



dhW(x,t) 
Jt 



2tt 

+ 8 ri dy 
3 is ?/ 



log(l 



x , 



/i (2) (x,t) 



1-2/ 



h^(-,t)-yh^\x,t) )-hW(-,t) 

y y 



(30) 



As indicated above, the structure functions must vanish at the boundary x — 1 in order 
to cancel the divergence of the logarithm in eq (^) and thus for the GLAP procedure to 
be applicable. This makes the projection of the rest frame structure functions mandatory. 
The variation of the structure functions for finite intervals in t is straightforwardly ob- 
tained by iteration of these equations, i.e. as a solution to the differential equation (j30"l). 
As discussed previously the initial value for integrating the differential equation is given 
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by the scale Ql at which the model is defined. It should be emphasized that this scale es- 
sentially is a new parameter of the model. For a given constituent quark mass m we adjust 
Ql to maximize the agreement of the predictions with the experimental data on previously 
[ |T7f calculated unpolarized structure functions for electron-nucleon DIS: — F™ ■ For 
the constituent quark mass m = 400MeV we have obtained Ql ~ 0.4GeV 2 . Note that this 
value of Qq is indeed (as it should) smaller than the ultraviolet cut-off of the underlying 
NJL soliton model as A 2 m 0.56GeV 2 . The latter quantity indicates the range of validity 
of the model. In figure |5|a we compare the un-evolved, projected, proton structure func- 
tion h p T (x, Ql) with the one evolved from Ql = 0.4GeV 2 to Q 2 = 4.0GeV 2 . As expected 
the evolution pronounces the structure function at low x. 

This change towards small x is a generic feature of the projection and evolution pro- 
cess and presumably not very sensitive to the prescription applied here. In particular, 



choosing a projection technique |33j alternative to ( Pq) may easily be compensated by an 
appropriate variation of the scale Qq. In figure |]b the same calculation for (x,Q 2 ) is 
presented. 

In the evolution of the twist-2 pieces we have restricted ourselves to the leading order 
in a s because for the twist-3 piece of Hl, the necessary ingredients are not known in 
next-to-leading order. Even the leading order evolution is only known in the large Nc 
limit. It should be noted that such an approach seems particularly suited for soliton 
models which equally utilize large Nc arguments. As pointed out by Balitskii et al. 



34f the admixture of independent quark and quark-gluon operators contributing to the 



twist-3 portion Hl(x,Q ) grows with n where n refers to the n moment, A4 n h^iQ 
of hi{x,Q 2 ). However, much like the case with the spin-polarized structure function, 
g 2 (x,Q 2 ) |35| in the Nc — ► oo limit the quark operators of twist-3 decouple from the 
quark-gluon operators of the same twist. Then the anomalous dimensions 7„ which 
govern the logarithmic Q 2 dependence of M. n can be computed. Once the 7 n 's are known 
an evolution kernel can be constructed that "propagates" the the twist-3 part h(x, Q 2 ) 
in momentum 

M^Q 2 ) = j 1 ^ L b{x,y ] Q 2 ,Ql)h L {y,Ql) ■ (31) 
Jx y 

We relegate the detailed discussion of the kernel b(x,y;Q 2 ,Ql), which is obtained by 
inverting the Q 2 dependence of M. n , to appendix C. In figure |]a we show the evolution 
of h L (x). Again we used Ql = 0.4GeV 2 and Q 2 = 4.0GeV 2 . 

As discussed in ref the merit of this approach is that to leading order in Nc 
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the knowledge of hi{x,Q 2 ) at one scale is sufficient to predict it at any arbitrary scale, 
which is not the case at finite Nc- 5 Thus Iil{x,Q 2 ) obeys a generalized GLAP evolution 
equation. This finally enables us (in much the same manner as was the case for g?,{x, Q 2 ) in 
|p.8|| ) to compute the longitudinal chiral odd structure function hi(x, Q 2 ) by combining the 
separately evolved twist-2 and twist-3 components together. The result for Ql = 0.4GeV 2 
and Q 2 = 4.0GeV 2 is shown in figure ^|b. We recall that the only ingredients have been 
the leading twist pieces of the chiral odd structure functions at the model scale Qo-f\ 



5 Discussion of the Numerical Results 



In this section we discuss the results of the chiral-odd structure functions calculated from 
eqs (|59l) -(|62"D for constituent quark masses m = 400MeV and m = 450MeV. In figure [l] we 
have shown the up and down quark contributions to the transverse chiral odd structure 
function of the proton. Figure [| displays them boosted to the IMF. We observe that 
these structure functions are always smaller (in magnitude) than the twist-2 polarized 
structure function g\ with the same flavor content. This relation is also known from the 
bag model 0] . Similar to the confinement model calculation of Barone et al. we find 
that hj . (x) is negative at small x. In contrast to g[ d \x), however, it might change sign 
although the positive contribution appears to be small and diminishing with increasing 
constituent quark mass. 

As already indicated in the introduction the DIS processes which are sensitive to 
these distributions will provide access to the charge weighted combinations thereof. We 
will hence concentrate on this flavor content. In any event, as we will be discussing 
both, the proton and the neutron chiral odd distributions, other flavor combinations can 
straightforwardly be extracted by disentangling the isoscalar and isovector pieces in eq 
fllTf). In connection with the chiral-odd transverse nucleon structure function we also 
calculate its zeroth moment which is referred to as the isoscalar and isovector nucleon 
tensor charges ||], 

TUQ 2 ) = y [ [dx h p T (x, Q 2 ) + h n T {x, Q 2 )] (32) 



s As noted in pi| , next to leading order corrections are estimated to go like O (l/N^ x ln(n)/n) at 
large n. 

6 A feature of h^(x) compared with (x) is that as hi, (x) does not mix with gluon distributions owing 
to its chiral-odd nature and its Q 2 evolution is given by (pa), (|6^ ) even for the flavor singlet piece. 
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T V T {Q 2 ) = 6 



dxK p T (x,Q 2 ) -h%(x,Q' 



(33) 



at both the low scale, = 0.4GeV 2 and a scale commensurate with experiment, Q 2 = 
4GeV 2 . Of course, for the neutron we have to reverse the signs of the isovector pieces 
in eq (0). I* 1 ec L s (H3) anc ^ dUD the normalization factors are due to the separation into 
isosinglet and isovector contributions, cf. eq (]I9|) . Note that due to Jq dzP^ g (z) ^ 
the tensor charge is not protected against logarithmic corrections. Our results for the 
valence quark approximation are summarized in Table 1. For completeness we also add 
the vacuum contribution to the tensor charges at the model scale Qq. Their analytic 
expressions are given in appendix D. Obviously this vacuum contribution is negligibly 
small. This is a strong justification of the valence quark approximation to the chiral odd 
structure functions. A further justification comes from a recent study of the Gottfried 



Table 1: Nucleon tensor charges calculated from eqs (p2|) and (p3|) as a function of the 
constituent quark mass m in the NJL chiral-soliton model. The momentum scales are 
<5o = 0.4GeV 2 and Q 2 = 4.0GeV 2 . The numbers in parenthesis in the respective upper 
rows include the negligible contribution from the polarized quark vacuum. We compare 
with results from the Lattice P6]| , QCD sum rules JTjJ], the constituent quark model 
with Goldstone boson effects [Q and a quark soliton model calculation |3£| including 
multiplicative 1/Nc corrections violating PCAC in the similar case of the axial vector 
current f40R . Finally the predictions from the confinement model of ref || with the 
associated momentum scales (in GeV 2 ) are shown. 



m (MeV) 


350 400 450 


Lat. SR CQ QS 


Q 2 CM 


rf(Q 2 ) 


0.80 (0.82) 0.72 (0.76) 0.67 (0.72) 
0.73 0.65 0.61 


0.61 0.61 1.31 0.69 
no scale attributed 


0.16 0.90 
25.0 0.72 




0.88 (0.89) 0.86 (0.87) 0.86 (0.85) 
0.80 0.78 0.77 


1.07 1.37 1.07 1.45 
no scale attributed 


0.16 1.53 
25.0 1.22 



sum rule within the same model 11411. Also in that case the contribution of the distorted 



quark vacuum to the relevant structure function turned out to be negligibly small. 

Besides justifying the valence quark approximation for the chiral odd distributions 
table [l] contains the comparison to other model calculations of the nucleon tensor charges. 
We note that in obtaining the isovector tensor charge Tj, we have omitted contributions 
which are suppressed by 1/Nq {cf. appendix D). These contributions arise when one 
adopts a non-symmetric ordering of the operators in the space of the collective operators 
[p9| . The main reason for taking the symmetric ordering is that in the case of the isovector 
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axial charge, g^, any non-symmetric ordering of the collective operators leads to a sizable 
violation of PCAC unless the meson profile is not modified p0 |. These multiplicative 
1/Nc corrections [[E| may be the reason why our predictions for are somewhat lower 
than those of other models. In the case of the flavor singlet component, which does not 
have such corrections, our results compare nicely with other model calculations except for 
the constituent quark model of ref |38|| . 

In figure [3]we display the transverse chiral odd proton \)? T (x, Qq) and neutron (x, Ql) 
structure functions at the low momentum scale Q$, while in figure |] we do the same for 
the corresponding chiral odd longitudinal structure functions h p L (x,Ql) and h\ (x, Q%). 



0.45 



m = 400 MeV 
m = 450 MeV 



0.15 



0.05 



m = 400 MeV 
m = 450 MeV 




0.00 



3.05 



Figure 3: The valence quark approximation of the chiral-odd nucleon structure functions 
as a function of Bjorken-x. Left panel: hP T (x, Ql) for constituent quark masses m = 
400MeV (solid line) and m = 450MeV (long-dashed line). Right panel: (x,Ql). 

We observe that the structure functions h^(x, Ql) and (x, Ql) are reasonably localized 
in the interval < x < 1. In particular, this is the case for the chiral odd structure 
functions of the neutron. Nevertheless a projection as in eq (EBp is required to implement 
Lorentz covariance. In addition the computed structure functions exhibit a pronounced 
maximum at x ~ 0.3 which is smeared out when the constituent quark mass m increases. 
This can be understood as follows: In our chiral soliton model the constituent mass serves 
as a coupling constant of the quarks to the chiral field (see eqs (|]) and (§)). The valence 
quark level becomes more strongly bound as the constituent quark mass increases. Hence 
the lower components of the valence quark wave-function increase with m and relativistic 
effects become more important. This effect results in the above mentioned broadening of 
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Figure 4: The valence quark approximation of the chiral-odd nucleon structure functions 
as a function of Bjorken-x. Left panel: h p L (x, Qq) for constituent quark masses m = 
400MeV (solid line) and m = 450MeV (long-dashed line). Right panel: h n L (x, Qq). 



the maximum. 

As discussed above a sensible comparison with (eventually available) data requires 
either to evolve the model results upward according to the QCD renormalization group 
equations or to compare the model results with a low momentum scale parameterization 
of the leading twist pieces of the structure functions. The latter requires the knowledge of 
the structure functions at some scale in the whole interval x G [0, 1[. At present no such 
data are available for the chiral odd structure functions Ht{x) and Jil(x). Therefore and 
in anticipation of results from RHIC and or HERMES we apply leading order evolution 
procedures to evolve the structure function from the model scale, Qq = 0.4GeV 2 to 
Q 2 = 4GeV 2 . In Figs. [5|a and |5|b we display the results of the two step process of 
projection and evolution for the twist-2 transverse structure function, hj>(x,Q 2 ) and 
h p ^ 2 \x,Q 2 ), respectively for a constituent quark mass of m = 400MeV. In figure || we 
present the evolution of h p L (x) along with its decomposition into terms of the leading 
twist-2 contribution, 2x dyhq-(y Q 2 ) / 'y 2 \ and the remaining twist-3 piece, li^x^Q 2 ). 
As in the case of the polarized structure function, g2{x, Q 2 ), the non-trivial twist-3 piece 
arises as a result of the binding of the constituent quarks through the pion fields acting as 
effective non-perturbative gluonic modes. The twist-3 contribution is evolved according 
to the large Nc scheme []33], |3| outlined in the preceding section (and in Appendix 



C). Similarly in Figs. [7] and [8] we display the projection and evolution procedure to the 
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Figure 5: Left panel: The evolution of K P T (x, Q 2 ) from Ql = 0.4GeV 2 (solid line) to 
Q 2 = 4GeV 2 (long-dashed line) for the constituent quark mass m = 400MeV. Right 
panel: The evolution of the twist-2 contribution to the longitudinal chiral odd structure 
function, h p L {2) (x, Q 2 ) from Q 2 = OAGeV 2 (solid line) to Q 2 = 4GeV 2 (long-dashed line) 
for m = 400MeV. 

twist-2 and 3 contribution to the neutron structure functions, fi^ 2 \x, Q 2 ) and 7^(x, Q 2 ), 
respectively. 

Besides the absolute magnitudes, the major difference between the chiral odd structure 
functions of the proton and the neutron is that the latter drop to zero at a lower value of 
x. As can be observed from figure |]this is inherited from the model chiral odd structure 
function at the low momentum scale and can be linked to the smallness of the down quark 
component of tiT, cf. figure 0. Apparently the projection and evolution program does not 
alter this picture. 

We would also like to compare our results from the NJL chiral soliton model to those 
obtained in other approaches. A MIT bag model calculation of the isovector contribution 
6(fij, — hJf) has been presented in ref 0]. In shape {e.g. position of the maximum) that 
result is quite similar to ours. However, the absolute value is a bit larger in the MIT bag 
model. This reflects the fact that in the MIT bag model the isovector combinations of the 
axial and tensor charges turn out to be bigger than in the present model. Additionally, the 
QCD evolution of the MIT bag model prediction for Ht has been studied in ref utilizing 
the Peierls-Yoccoz projection as in ref |21|| . In that case the maximum at x ~ 0.5 gets 
shifted to a value as low as x = 0.2. Also the structure function becomes rather broad at 
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h^(x,Q 2 ) 
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Figure 6: Left panel (|a): The evolution of the twist-3 contribution to the longitudi- 
nal chiral odd structure function, h P L (x,Q 2 ) along with the corresponding twist-2 piece, 



h p L {2) (x,Q 2 ). Right panel (|b): The evolution of h\ (x,Q 2 ) = h p L (2) (x,Q 2 ) + V L {x,Q 



P(2) 



from Q 2 = 0.4GeV 2 (solid line) to Q 2 
quark mass m = 400MeV. 



4GeV (long-dashed line) for the constituent 



the large scale. The fact that in that calculation the evolution effects are more pronounced 
than in the present approach is caused by the significantly lower scale (//bag = 0.08GeV 2 ) 
used in ref ||. On the other hand our results are quite different to those obtained in the 
QCD sum rule approach of ref [ft3f| . The sum rule approach essentially predicts h T to be 
constant in the interval 0.3 < x < 0.8. For small values of x the authors of ref |43| assume 
a Regge behavior. In the (covariant) constituent quark model of Suzuki and Shijetamin 



[p8| a result similar to ours is obtained when effects attributed to Goldstone bosons are 
included. Otherwise the maximum of their distribution is about 50% larger than in our 
calculation. These authors also observe that in magnitude the down quark component 
is significantly smaller than the up quark piece. The chiral chromo-dielectric model of 
Barone et al. [[J predicts a similar shape for hr but their distribution is larger than 
the one in NJL chiral soliton model. This is also reflected by the sizable value for the 
isovector tensor charge T^(Q 2 = 25GeV 2 ) = 1.22 and T^(Q 2 = 0.16GeV 2 ) = 1.53 in that 
approach. In the Isgur Karl model (which has hi — Qi) the maximal value of h\ is only 
about have as big as in our calculation, cf. figure 1 in ref 0. 

For completeness we also demonstrate in figs U and ITUI that at the low model scale 



Ql = 0.4GeV 2 Soffer's inequality |44]] is satisfied. This inequality relates the nucleon chiral 
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Figure 7: Left panel:The evolution of hSj. (x, Q 2 ) from Qq = 0.4GeV 2 (solid line) to 
Q 2 = 4GeV 2 (long-dashed line) for the constituent quark mass m = 400MeV. Right 
panel: The evolution of the twist-2 contribution to the longitudinal chiral odd structure 
function, hf 2) (x, Q 2 ) from Q 2 = OAGeV 2 (solid line) to Q 2 = 4GeV 2 (long-dashed line) 
for m = 400MeV. 



odd distribution functions to both the unpolarized f\ (x, Qq) an d polarized g[ q> (x,Ql 



structure functions 



fl q \x 1 Q 2 )+g?\x,Q 2 )>2 \hi°>(x,Q 



(9), 



(34) 



Here the superscript refers to the flavor combination which projects onto up and down 
quark quantum numbers (q = u,d). Note again, that this projection refers to the con- 
stituent quarks which contain some non-perturbative gluonic distributions. In figure TG 
we display the down quark component of the inequality fl3"4"|). As this component of Ht 
is almost negligible the inequality is satisfied by the unpolarized structure function 
being larger in magnitude than the polarized one g[ . For the constituent quark mass 
m = 450MeV we also find that Soffer's inequality is satisfied at the model scale Qo- The 
only remarkable difference to m — 400MeV is that the down quark component of hr? 
does not possess nodes. This can already be inferred from figures [l] and [| 

A more thorough model study of Soffer's inequality for scales other than Qq which 
also contains the next-to-leading order contributions in the evolution program is subject 
to further investigations [45]. The next-to-leading order calculation for in the 
Isgur-Karl model indicates that its scale dependence is slightly mitigated by the inclusion 
of next-to-leading order contributions. 
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Figure 8: Left panel: The evolution of the twist-3 contribution to the longitudinal 
chiral odd structure function, h L (x,Q 2 ) along with the corresponding twist-2 piece, 
hl {2) (x,Q 2 ). Right panel: The evolution of h n L (x,Q 2 ) = h n L {2) (x,Q 2 ) + 7^(x,Q 2 ) from 
Ql = 0.4GeV 2 (solid line) to Q 2 = 4GeV 2 (long-dashed line) for the constituent quark 
mass m = 400MeV. 

6 Conclusions 

In this paper we have presented the NJL chiral soliton model calculation of the leading 
twist parts of the transverse and longitudinal chiral odd structure (distribution) functions 
of the nucleon. Data on these distribution functions should eventually be available from 
DIS experiments in the fragmentation regions (in conjunction with fragmentation func- 
tions) or be extracted from Drell-Yan experiments. These structure functions serve to 
complete our picture of the spin distributions of the nucleon. The most important feature 
of the present quark based model is that it is chirally invariant and that this symmetry 
is dynamically broken. After bosonization the NJL model becomes an effective meson 
theory in which baryons emerge as self-consistent soliton solutions exactly the way as 
expected from large N c considerations in QCD. Chiral soliton models are particularly 
interesting in the context of the nucleon spin structure as these models nicely explain the 
small contribution of the quarks the total nucleon spin. 

In the NJL chiral soliton model there are two contributions to nucleon properties. 
First, there is the contribution of the distinct valence quark level. This is the lowest level 
in the quark spectrum and bound in the background of the chiral soliton. Second, there 
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Soffer Inequality 



Soffer Inequality 




X X 

Figure 9: Left panel (§a): The Soffer inequality for the chiral even combination 
fi (x, Qq) + g[ u \x, Qq) (solid line) of the effective up-quark distributions and the chiral 
odd structure function 2 h^\x,Ql) (long-dashed line) for a constituent quark mass of 
m = 400MeV, calculated in the nucleon rest frame (RF). Right panel: (§b) Same as figure 
§a calculated in the infinite momentum frame (IMF). The transformation prescription is 
given in eq ([2~6|). 

is the part which is associated to the polarization (by the soliton) of the vacuum. For 
many static nucleon properties the latter contribution is quite small, in particular for 
those which are related to the axial (spin) properties of the nucleon. This is a strong 
indication that the vacuum contribution to the chiral odd structure functions is negligible 
as well. Hence it seems more important to include substantial 1/Nc corrections to the 
valence quark contribution. These corrections come about when projecting the soliton 
onto states with good spin and isospin, i.e. proton and neutron. Inclusion of these 1/Nq 
corrections together with a consistently regularized treatment of the vacuum polarization 
is technically rather involved and beyond the scope of the present paper. The numerical 
results for the tensor charge indicate that the vacuum contributions are in fact negligibly 
small. 

When the model structure functions are computed one immediately recognizes that 
they have improper support due to the breaking of translational invariance by the back- 
ground soliton, i.e. the structure functions do not exactly vanish for x > 1. This can be 
cured by Lorentz boosting to the infinite momentum frame which is particularly suited 
for DIS processes. Although the un-boosted structure functions are negligibly small at 
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Figure 10: Same as figure [5] for the down quark combination. 



x > 1 the transformation to this frame is essential and has sizable effects on the struc- 
ture functions at moderate x. However, the most important issue when comparing the 
model predictions to (not yet available) experimental data is the observation that the 
model represents QCD at a low momentum scale Qq. A priori this scale represents an 
additional parameter to the model calculation which, for consistency, has to be smaller 
than the ultraviolet cut-off of the model A 2 = 0.56GeV 2 . For the model under considera- 
tion we previously fixed Ql when studying the unpolarized structure functions and found 
Ql = 0.4GeV 2 . The important logarithmic corrections to the model structure functions 
are then obtained within a generalized GLAP evolution program. In this context we have 
restricted ourselves to a leading order (in «qcd) calculation because the anomalous di- 
mensions, which govern the QCD evolution, for the twist-3 piece of the longitudinal part 
of the chiral odd structure are only known to that order. As the full evolution to the 
longitudinal structure function involves both twist-2 and twist-3 pieces this restriction 
is consistent. We have seen that the QCD evolution of the chiral odd structure function 
leads to sizable enhancements at low x, i.e. in the region 0.01 < x < 0.10. In this respect 
the present situation is similar to that for the polarized structure functions. A difference 
to the polarized structure function is that the lowest moment is not protected against 
logarithmic corrections, even at leading order in ttQCD- For the nucleon tensor charge 
we thus find a reduction of about 10% upon evolution to Q 2 = 4.0GeV 2 . We have also 
compared the neutron and proton chiral odd structure functions. This has been achieved 
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by the inclusion of the 1/Nc cranking corrections. In absolute value the proton struc- 
ture functions are about twice as large as those of the neutron. Furthermore the neutron 
structure functions drop to zero at a lower value of x. These two effects can be linked 
to the down quark component of the transverse nucleon chiral odd distribution functions 
being significantly smaller than the component with up-quark quantum numbers. We 
have also observed that neither of these features is effected by the evolution program. 
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A Bilocal Light Cone Distributions 

In this appendix we outline the steps giving rise to the starting point of our calculation, 
eqs ( [RD and (|T^) in section 3. 

It is well known that the chiral odd spin-dependent structure functions hxix) and 
hi(x) do not contribute to the hadronic tensor in DIS. hr{x) was first studied by Ral- 
ston and Soper in the context of polarized Drell-Yan processes while hi{x) was more 
recently detailed by Jaffe and Ji Q. In the latter study a general Lorentz decomposi- 
tion of invariant matrix elements of the characteristic bilocal operators in hard processes, 
\I / (0)r At \E'(A?2) was performed.^] Adopting light-cone variables reveals that up to twist— 3 
there are six invariant structure functions that characterize the nucleon. The leading twist 
(two and three) spin-dependent contributions emerging from the Lorentz decomposition 
of these hard processes are given by||, 

/ — e* Xx (P, (0)icv7 5 * (An) \P, S) = 2 {h T (x) (S XfJ p v - S ±vPli ) / M 

+ h L (x)M (p^p u - pup,,,) S • n} , (35) 

7 See ref |43| for the definition of chiral odd structure functions in the language of a hadronic tensor. 
Actually such a definition is sufficient to carry over the QCD definition of the chiral odd structure 
functions to the NJL model because formally the NJL model currents are identically to those in QCD. 
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/ — e iXx (p, S|*(0)z7 M 75* (An) \p, S) = 2 {g L (x)p,S ■ n + g T (x)S ± ,} (36) 

where p M and n M define a light-like coordinate system, i.e. p ■ n = 1, p 2 = n 2 = and 
IPS') denotes the nucleon state with four momentum P and spin S. In the system where 
the nucleon is moving along the z direction one conveniently defines the four vectors 

P" = T^ 1 ' ' ' 1 ) 

n» = -^P(l,0,0,-1) . (37) 

In this system the nucleon momentum is given by P^ = p^ + Mn fl /2 and spin is 
decomposed as = (S ■ n)p^ + (S ■ p)^ + Sj_. Finally, V — > oo corresponds to the 
infinite momentum frame (IMF) and V = Mj v2 corresponds to the nucleon rest frame 
(RF). Utilizing the convenient projection properties of the light-like vectors the defining 
equation (|3~5D may be inverted. One obtains for the chiral odd structure functions 

1 r d\ 



h T (x) = -^J^e^(PS ± \^ + (0) 1± ^ + (Xn)\P,S ± ) (3* 



and 



h L (x) = ^J^(PS z \yl(0) lo ^ + (\n) 

-^ t +(0)7o75*-(A^)|P,S 2 ) . (39) 

The quark bilocals describe the propagation of the intermediate constituent quark which 
is struck by the external source. The forward propagation is described by x > while 
negative x parameterize an intermediate quark which moves backward. In what follows we 
will only consider positive x in conjunction with the contribution associated with the for- 
ward propagating quark h^ + \x). The backward contribution can easily be obtained from 
h^ + '(— x). Finally, noting the change of variables from light-like coordinates (77, A, to 
light-cone coordinates (£ + ,£~,£j_) where in particular, 

t = r)V and r = ^ (40) 

yields the chiral odd quark transverse distributions, 

hT~\x) = j dOsxpHf - ^) 

x (S ± \^ + (Ol±^ + (0)\S±) c+ _ t _ n , (41) 
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and the longitudinal contribution, 

x(^|* t + (07o75*-(0) - ¥_(£) lo j 5 Q 2 y + (0)\Sz)t+=^=o • (42) 
These equations represent the starting point of section 3. 

B Chiral Odd Structure Functions in the NJL Soli- 
ton Model 

In this appendix we derive and summarize the explicit expressions for the chiral odd 
structure functions, eqs The first step is to construct the eigenfunctions of the single 
particle Dirac Hamiltonian (|8|) in coordinate space. The hedgehog ansatz ([?]) connects 
coordinate space with isospace and these eigenfunctions are also eigenstates of the grand 
spin operator 

G = J+ i = ' + f + i («) 

which is the sum of the total spin J and the isospin r/2. The spin itself is decomposed 
into orbital angular momentum / and intrinsic spin er/2. Denoting by M the grand spin 
projection quantum number the tensor spherical harmonics which are associated with the 
grand spin may be written as yf- M (r). Note that these tensor spherical harmonics are 
two-component spinors in both spin and isospin spaces. Given a profile function G(r) 
the numerical diagonalization of the Dirac Hamiltonian (§) yields the radial functions 
fiu G '' + ' 1 ' ) ( r )? / / i G ' + ' 1 ' ) ( r ); e tc- in the decomposition (c/. ref p6|) 



2 

The second superscript (±) denotes the intrinsic parity, which also is a conserved quantum 
number.0 Note that for the G = channel, which contains the mean-field contribution 
to the valence quark wave-function in eq (O) 



wAr)y°A (f) 



3 The total parity is given by the product of the intrinsic parity and (— ) c 
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only the components with j = + 1/2 are allowed. In addition to this mean-field piece 
(f46|) the complete valence quark wave-function ( |i~3D also contains the cranking correction, 
which dwells in the channel with G = 1 and negative intrinsic parity. 

The discretization (//) is accomplished by choosing suitable boundary conditions at a 
radial distance which is large compared to the soliton extension [[5], H7J. This calculation 
yields the energy eigenvalues e M , which enter the energy functional @. The soliton con- 
figuration is finally determined by self-consistently minimizing this energy functional. In 
ref ES[ the numerical procedure is described in detail. 

We continue by making explicit the Fourier transform of eq (|13|), 

~ f d^x ~ ~ 

(p) = J -^Vv (x) exp (ip ■ x) = * v (p) + „ (p) . (47) 

The leading order in iV<7 valence quark contribution is just the Fourier transform of (|46| ) 

*" (p) = ! ( /, M ^ ») J (48) 

and the cranking correction involves the Fourier transform of spinor with G = 1 and 
negative intrinsic parity 

%(P) = - ( fir, w ^ _ /f w ^ ^ j • (49) 

Here yf- M (ft) are the Fourier transforms of the tensor spherical harmonics associated 
with the grand spin operator fl4"3|). The Fourier transform for the radial functions in eqs 
(H) and (H) is defined by 

rR 

4>Jp) = dr r 2 ji(pr)(j)^(r) . (50) 



o 

Here the index I of the spherical Bessel function denotes the orbital angular momentum 
of the associated tensor spherical harmonic. We have suppressed the grand spin index 
on the transforms of the radial wave functions for convenience. For purposes of notation 
we have also introduced the quantity in eq (57T) which parameterizes the cranking 
corrections in eq ( |TB| ) 

(/i|T '" |v) - o?Q, fe^+^-^^-^-^o}^ (51) 

Q, fe (Jx + ih) - ^ (Ji - ih) ~ S Mfl J \ *g Mi i (52) 
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Q, = T{ l _ \ i dr r * & (r) gf ) (r) + / v (r) /f (r)} . (53) 



where 



a 2 (e, 

In this definition we have included the total moment of inertia a 2 . In the proper-time 
regularization of the NJL chiral soliton model a 2 is given by |L9| 



^ E lfrWf, (54) 



a 



2 AT C 



2 ^^(A)(/,|r»^|r» , (55) 
a 2 = ^(l + sgn(e val ))a 2 + a s 2 . (56) 



The regulator function in the vacuum contribution reads 



A e -(WA) 2 _ e -(^M) 2 sgn(e„)erfc (|x|) ~ sgn(e M )erfc ( 



«A) = -t=- -5—3 — ^/ — • (57) 



The moment of inertia enters via the quantization description for the collective coordinates 
Q — > a 2 ./ with ,7 being the nucleon spin operator. In this quantization prescription we 
had previously restricted the moment of inertia to its valence quark contribution, a 2 , to 
ensure that the Adler sum rule for the unpolarized structure functions is maintained in 
the valence quark approximation [17]. For small or moderate constituent quark masses 
the valence contribution to the moment of inertia is about 80% or more |E7| . This is one 



of the reasons to believe that the valence quark approximation to structure functions is 
sensible. In the case of the chiral odd structure functions (as for the polarized ones) the 
valence quark approximation appears to be even better. As we see from table [l] the lowest 
moments of these structure functions are saturated to about 95% by the valence quark 
contribution. Hence it is reasonable to assume that the vacuum contribution to these 
structure functions in negligibly small. As a consequence the valence quark approximation 
with the total moment of inertia substituted into the quantization rule will provide a very 
reliable estimate of the chiral odd structure functions. 

Together with (N\Dij\N) = — (4/3)IjJj [|13| the nucleon matrix elements may now 
easily be computed. Here J denotes the nucleon isospin. Whenever products of collective 
coordinates and operators appear which do not commute after canonical quantization we 
adopt the symmetric ordering. This is consistent with fundamental requirements such as 
PCAC. Defining finally the following combinations 

/ W (P) = Qj ( ; ] (P) and (p) = Q,~gf (p) , (58) 
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for i = 1,2. The isoscalar(vector) contributions to the chiral odd structure functions fl24|) 
read 



N, 



5M 



N 



367T JM N \x T \ 



pdp 



L , , , . 3cos 2 (^)-l 1. (2) 
x j# v (p)g W (p) ^ (p) 



± (p) (p) + /v (p) (P) 

=F (/v (P) ^ (P) + (P) / (2) (P) 



cos(0±) 



cos 



-/v (P) (P) 



cos 2 (^)-3 - / N ?(a) ^ 008(^2' 



4v/2 



" /v (P) / W (P) 



(59) 



iv. 



M; 



N 



C 



pdp 



367T JMiv|x T | 

x {g v (pf ± 2^ v (p) / v (p) cos(0±) + / v (p) 2 (cos 2 (^)) } 



(60) 



N, 



5M 



TV 



367T Jm n \x=f\ 



pdp 



X {±^v (P) ^ (P) 3COS2 4 ( ^ 1 ± ~<7v (P) ^ (P) 



=F/v (P) (P) 



; 1 ^ ' ^ 1 + C ° S2( ^ ) T /v (P) / (2) (P) 



2cos 2 (^) - 1 



2v/2 



(61) 



A/7" poo 

h%t(x, fi 2 ) = -Nc^ J Pdp {t~ 9v (P) 2 ± A (P) 2 (2cos 2 (^) - l) } . (62) 



which we evaluate numerically. Note that the angle 9^ is related to the integration variable 
p via 



cos^ = -\MpfX ± e v | 
^ p 



(63) 



In ref |TIJ the contribution to structure function hx from effective quark distributions^ 
was calculated omitting the cranking corrections and adopting an external (non-self- 
consistent) meson profile. 



9 In this work it is important to note that the quark distributions refer to constituent quarks, m q rs 
400MeV, it is thus misleading to compare them with the data of parton distributions from either Drcll- 
Yan or DIS processes. 
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C Evolution of 1il(x, Q a 



In this appendix we outline our technique to evolve the low scale model prediction for the 
twist-3 piece 1il(x, Q%) to the larger scale Q 2 . This utilizes the method described in refs 
P, ffUJ based on the results of ref [EDI. The Q 2 evolution of the moments 



M n h L (Q* 



is given by 



M n h L (Q 2 



dxx n hi J (x 1 Q 2 



'a(Q 2 )V" /bo r 

W ] » M n h L (Ql) 



(64) 



AQl) 



(65) 



Here 60 — (UNc — 2n/) /3 is coefficient of the leading term in the QCD-beta function. 
Also, Nc and n/ are the number of colors and flavors respectively. Within the Nc — > 00 



approximation the anomalous dimensions are [34 



7 £ = 2N C [S n + 7e - - + 



4 2n+l 



(66) 



with 5 n = X)jLi(l/j) — 7e where 7^ = 0.577. . . is the Euler constant which has been 
introduced for later convenience. 

In order to find the QCD-evolution of the structure functions one needs to invert the 
Mellin-transform fl6"4"|). This can be achieved by noting that the Bernstein polynomial 



(67) 



has the property, 

lim b (N ' k) {y) = 5{x-y) 

N,k^>oo 

for < x,y < 1. This enables one to express the structure function via its moments 

(N + iy. N ^ k (-iy 



(6* 



h L (x,Q') 



lim 

k/N-+x 



E 



% l\(N-k-l)\Jo 



dyy k+l h L (y,Q' 



(69) 



which depend on Q 2 as indicated in eq 



h L (x,Q 2 



(iV+i)!^ 
lim > 

Af,fc^oo A;! f-' 



-11 



l\(N — k — 1)1 



i>V 6 ° ( dyy k+l h L (y,Ql) . (70) 
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Here L = a(Q 2 )/a(Ql) denotes the ratio of the running coupling constants in QCD. 
Unfortunately, the rapid oscillations in the summation over I in (|70|) due to the factor 
(—1)' preclude numerical summation of (|70P . Yet, observing that the expression L 7 ™^ 
may be expanded as 

Cj(L) 
. ;0 (n+py-rW 



L^ b ° = a{L) 



(71) 



where a(L), and r(L) are constants determined from the asymptotic form (n — > oo) of eq 



r{L) = 2N c ln(L)/b and a(L) = exp 



r(L) ( 7e - ^ 



(72) 



one can perform the sum to any desired accuracy. It should be noted that p remains 
undetermined. It may be varied to control the convergence of the series (ffl|) . To determine 



the expansion coefficients Ci(L) we rearrange eq ( 7lD to a Fourier expansion, 



^1 — zp) r exp 



1 



2n 2n + 2 



E 

k=l 



2k 



2kn 2k 



(73) 



i=0 



Here z = l/(p + n)^n=l/z — p and Ci(L) = Ci(r). Furthermore we have utilized the 
asymptotic expansion of 

1 ^ / B 2k 



S n = ln(n) + 



2n 



E 

fc=i 



\2kn 2k 



(74) 



where the ^fc's are the Bernoulli numbers. Performing a Taylor series to eighth order in 
z yields the following values for the expansion coefficients Ci(L) = Cj(r(L)) (for p — 2) 



C (r 
C 3 (r 
C 6 (r 
C 7 (r 
C*(r 



1, 
1 

- r, 
2 ' 



Ci(r) = 0, 
C 4 (r) 



Cofr) = — r, 



12 



61 



25 



r H r' 

120 288 



— r H r z 

2 24 



125 97 2 125 , 

r H r H r , 

252 288 10368 ' 

1 37 2 25 , 
- r H r H r , 

2 80 576 



121 



354341 2 665 ., 



240 r + 604800 ' 6912 



r° + 



625 



497664 

which gives more than adequate convergence of the series. 



Finally we may write 

h L (x,Q 2 ) 



—b(x,y;Q 2 ,Ql)h L (y,Q 2 ), 
x y 



(75) 



(76) 
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where, 

OMQl) = a(L) g^pJj^ (77) 

is the evolution kernel used in eq fl3ip. It has been gained by using the additional relation 

l im (N + l)\ N ^ k y k+l _ 9(y-x) (xV- 1 ( yy+f- 1 

jv,™oo k\ frt l\(N - k - l)\ (k + I + pY+p T(i + p)y\y) \x) 

k/N—+x 

For the numerical results presented in sections 4 and 5 we have verified the stability of 
this evolution procedure by varying the undetermined parameter p in eq Q7T]) . 



D Tensor charges in the NJL chiral soliton model 

The conventional definition of the nucleon tensor charges reads 

(N\^a^\N) = rf ua^u and (N\^a^T 3 ^\N) = ua^r 3 u . (78) 

Here N again denotes the nucleon state. Note that both, the quark wave-function \l/ as 
well as the nucleon spinor u, are vectors in flavor space. Momentum labels have been 
omitted as the charges are defined at zero momentum transfer. Within the NJL chiral 
soliton model these charges can be extracted using standard techniques [Q: first, sources 



conjugated to the quark bilinears ^fa^ty and ^/a^ v r 3 ^ are added to the Lagrangian (|). 
Subsequently the bosonized and regularized action is expanded to linear order in both the 
sources and the angular velocities f2 (|IT|). The coefficients of the source terms then provide 
the charge operators in the space of the collective coordinates A, which are defined in eq 
(P~D|) . The corresponding matrix elements can be straightforwardly evaluated with the 
means provided in appendix B. Finally one obtains within the proper-time regularization 

r s N c (val|T»(H/?£ 3 T 3 |val) 

uj l v3i i tval t u 

Ac E^(A)(^k»(^S3r 3 |/i) (79) 



4a 2 

T V T = -^(l + S gn(e val ))(val|/3S 3 r 3 |val) 
o 



Nc 
6 



^( / u|/3S 3 r 3 |/x)sgn(e M ) erfc 



A 



SO) 
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Here \fi) denote the eigenstates of the static Dirac-Hamiltonian (]8|) and are the corre- 
sponding eigenvalues. Again |val) refers to the distinct valence quark level. The regulator 
function in the isoscalar piece (|79|) is identical to the one entering the moment of inertia, 
cf. eq (|57D. Those pieces containing the factor (1 + sgn(e va i)) are the valence contributions 
shown separately in table [L| 

As noted in section 5 we have omitted 1/Nc suppressed contributions to the isovector 
part which in the related case of the axial current violate PCAC 
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